Diffractive rj c and rjb productions by neutrinos via neutral currents 
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We report a first theoretical study for neutrino-induced diffractive productions of heavy pseu- 
doscalar mesons, rj c and T]b, off a nucleon. Based on factorization formalism for exclusive processes, 
we evaluate the forward diffractive production cross section in perturbative QCD in terms of the 
light-cone QQ wave functions (WFs) of T) c ,b mesons and the gluon distribution of the nucleon. The 
light-cone WFs of the r/ c (r/t) meson are constructed to satisfy the spin symmetry relations with 
those of the J /if) (T) meson. The diffractive rj c production is governed by the axial- vector coupling 
of the longitudinally polarized Z boson to QQ pair, and the resulting n c production rate is larger 
than the J/ip one by one order of magnitude. We also discuss the production of bottomonium n^,, 
which shows up for higher beam energy. 
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Exclusive diffractive leptoproductions of neutral vec- 
tor mesons provide unique insight into an interplay be- 
tween nonperturbative and perturbative effects in QCD 
[PJ. The diffractive processes are mediated by the ex- 
change of a Pomeron with the vacuum quantum num- 
bers, whose QCD description is directly related to the 
gluon distributions inside the nucleons for small Bjorken- 
x pi ||, |(| . The processes also allow us to probe the 
light-cone wave functions (WFs) of the vector mesons. 
Relating to the latter point, however, the applicability is 
apparently limited to probing the neutral vector mesons. 

In this Letter, we propose the exclusive diffractive pro- 
ductions of mesons in terms of the neutrino beam. The 
weak currents allow us to observe both neutral (M°) and 
charged mesons (M ± ) as iv + A -> v^j '/j^+N+M /M ± 
by Z/W boson exchange M , and these mesons can be not 
only vector but also other types of mesons including pseu- 
doscalar mesons. Thus, such processes may reveal struc- 
ture of various kinds of mesons, the coupling of the QCD 
Pomeron to quark-antiquark pair with various spin-flavor 
quantum numbers, and information on the CKM matrix 
elements. There already exist some experimental data for 
7T, p, Df, D* s §, §, D*+ @, and J/yj production |, 
but there are only a few theoretical calculations, e.g., for 
the J/ip production in a vector meson dominance model 
jll| and for D~ production with the generalized parton 
density p|. Our work gives a first QCD calculation for 
diffractive meson productions via the weak neutral cur- 
rent. Here, our interest will be directed to productions 
of heavy pseudoscalar mesons, especially rj c and r\b as in 
Fig. [I]. So far r\ c has been observed via the decays of J/ip 
or B mesons produced by pp and e + e _ reactions, while 
rjb has not been observed. The diffractive productions 
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via the weak neutral current will give a direct access to 
t] c as well as a new experimental method to identify rjb 
by e.g. measuring the two photon decay. 
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FIG. 1: A typical diagram for the exclusive diffractive t]q 
(Q — c, b) productions induced by neutrino (v) through the Z 
boson exchange. There are other diagrams by interchanging 
the vertices on the heavy-quark lines. 

We treat the rj c and rjb productions by generalizing the 
approach in the leading logarithmic order of perturbative 
QCD, which has been developed successfully for the vec- 
tor meson electroproductions |^|, |3|, ||, ^, |(| . We consider 
the near-forward diffractive productions Z* (q)+N(P) — > 
rj Q (q + A) + N'(P - A), where Q = c,b, and q,q + A,P 
and P— A denote the momenta of the virtual Z boson, tjq 
meson, initial and final nucleons, respectively. The total 
center-of-mass energy W — \/ (P + q) 2 is much larger 
than any other mass scales involved, i.e., W 2 ^> K 2 
with K 2 = Q 2 ,-t,m^,A| CD ,..., where Q 2 = -q 2 , 
t = A 2 , and tuq is the heavy-quark mass. We also sup- 
pose —t <C rriQ and Aq CD <§; to 2 -,. The heavy quark mass 
toq ensures that perturbative QCD can be applied even 
for Q 2 — to calculate the creation of the QQ pair as 
well as its time development before the nonperturbative 
effects to form a quarkonium state tjq set in. The crucial 
point is that at high W the scattering of the QQ pair on 
the nucleon occurs over a much shorter timescale than the 
Z* — > QQ fluctuation or the QQ — > tjq formation times 
(see Fig. [|). As a result, the production amplitudes obey 
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factorization in terms of the Z and j]q light-cone WFs. 
The QQ-N elastic scattering amplitude, sandwiched be- 
tween the Z and t)q WFs, is further factorized into the 
QQ-gluon hard scattering amplitude and the nucleon ma- 
trix element corresponding to the (unintegrated) gluon 
density distribution. This latter step in the factorization 
to get the gluon distribution can be carried out in the 
same manner as in the previous works J|, [| [|, q|. The 
participation of the "new players" Z and ijq in the for- 
mer step requires an extension of the previous works by 
introducing the corresponding light-cone WFs. 

First of all, we discuss the extension due to the par- 
ticipation of the Z boson. The ZQQ weak vertex of 
Fig. is given by (gw /2cos9 w )~fv( c v - ca7s), where 
(gw/2cos9w) 2 — V^GpM^ with Gf the Fermi constant 
and M z the Z mass. c v = 1/2- (4/3) sin 2 9 W , c A = 1/2 
for the c-quark and similarly for the 6-quark. As usual, 
we introduce the two light-like vectors q' and p' by 
the relations q = q' - (Q 2 /s)p', P = p' + {M 2 N /s)q', 
q' 2 = p' 2 = 0, s = 2q' ■ p' with Mm the nucleon mass and 
s = W 2 + Q 2 — Mjy, and the Sudakov decomposition of 
all momenta, e.g., k = aq 1 + f3p' + fcj_. We also introduce 
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FIG. 2: Diagrams representing the light-cone WFs: (a) for 
the virtual Z boson; (b) for the t]q meson. 

the polarization vectors e^' (£ = 0, ±1) of the virtual Z 

boson to satisfy ^(-l)^ 1 ^ £ ^ = -g^v+q^qu/M 2 z , 
which is the numerator of the propagator for the mas- 
sive vector boson. Because the q^q^/M^ term vanishes 
when contracted with the neutral current, we conve- 
niently choose as 



£ (o) = y + 

Q s 



5 < ±1 >=^ 1) = ^(0,l,±i,0),(l) 



for the longitudinal (£ = 0) and transverse (£ = ±1) po- 
larizations, respectively. Eq. ([!]) coincides with the polar- 
ization vectors for the virtual photon in Refs. J|, ||, |], |L 
For the J/-0 electroproductions, the wavy line of Figs, l] 
and ||(a) denotes a virtual photon. The corresponding 
photon light-cone WFs are derived as |, (§ |I| 



y/a(l-a)ux(k^v x >(q-k) 
a(l- a)Q 2 + k\+m 2 Q 



(2) 



with eq. (Q) for £ = 0, ±1. Here eq is the electric charge 
of the quark, and u\(k) (vx>(q — k)) denote the on-shell 
spinor for the (anti)quark with helicity A (A'). These 
WFs describe the probability amplitudes to find the pho- 
ton in a state with the QQ pair, where the quark carries 



the longitudinal momentum fraction a and the trans- 
verse momentum k± as k = aq' + /3p' + k^, k\ = —k\. 
Obviously, the light-cone WFs for the virtual Z bo- 
son can be obtained from eq. (|^) by the replacement 
e<27p -> (dw I '2 cos 9 W ) 7 p (cy - 0,475) as 



^\a,k 1 _) = - 



9w 



2 cos 9w 



\/a{l - a)u\(k)^ (c v - ca7s) vx>(q - k) 



a(l-a)Q 2 +k 2 



(3) 
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Next we proceed to the light-cone WFs for the t]q me- 
son, corresponding to Fig. ||(b) . Again, it is convenient to 
exploit the correspondence with the J/ip electroproduc- 
tions. The light-cone WFs for the heavy vector mesons 
V = J ftp, T have been discussed in many works, but are 
still controversial in the treatment of subleading effects 
like the Fermi motion corrections H j [l4|] , corrections to 
ensure the pure S-wave QQ state [[15 , corrections via the 
Melosh rotation |l6| , etc. Here we employ the light-cone 
WFs for the vector meson given by (see Fig. 0(b)) 



vx>(q - k) (zo) * n u x (k) (f>*(a,k ± ) 



M v 



,(4) 



where q = q + A, My , and = 0, ±1) are momen- 

tum, mass, and the polarization vector of the vector me- 
son with q 2 = M v , e^-q = 0, and e ^*-e (r °"> = -<Sw 
TZ = (1 + q'/M v )/2 denotes the projection operator, 
1Z 2 = TZ, to ensure the S'-wave QQ state in the heavy- 
quark limit (|, 14, M. (This projection operator coin- 
cides with that discussed in Rcf. (l5) up to the binding- 
energy effects of the quarkonia.) Note that eq. ([|) reduces 
to the vector-meson WFs of Ref. || by the replacement 
TZ — > 1. The physical interpretation of eq. (||) is simi- 
lar to the "perturbative" WFs (||) and (|^) , but the scalar 
function 4>{a, k±) contains nonperturbative dynamics be- 
tween Q and Q. Now, the light-cone WFs for the t\q me- 
son can be derived from eq. (|j) utilizing spin symmetry, 
which is exact in the heavy-quark limit. This symmetry 
relates the S'-wave states, t\q and the three spin states 
of the vector meson |l8[| . Namely, M nQ — My, and the 
pseudoscalar state is related to the vector state with lon- 
gitudinal polarization as \t\q) = 2Sq\V(zu = 0)), where 

Sq is the third component of the hermitean spin opera- 
tor Sq which acts on the spin of the heavy quark Q but 
does not act on Q. For the WFs we get (see eq. (|])) 



_ v\> (<? - k) 
vx< (q - 
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Here S 3 is a matrix repesentation of Sq as S 3 = 
7 5 ^°'/(2My), which is related to a spin matrix a 12 /2 = 
757°7 3 /2 in the meson rest frame by a Lorentz boost in 
the third direction Thus the r]Q is described by the 
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same nonperturbative WF <fr(a, k±) as the vector meson. 
Note that, due to the presence of TZ, the "QQtjq vertex" 
involves pseudovector as well as pseudoscalar coupling. 

Combining our Z and tjq WFs with the QQ-N elastic 
amplitude which was obtained in Ref . || , we get the total 
amplitude for the polarization £^ of the Z boson as 
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up to the terms suppressed for high IV and = — t — > 
0. Here a s (^j_) denotes the running coupling constant 
for V c colors. As usual, we have explicitly dealt with 
the imaginary part of the production amplitude, be- 
cause the small real part can be reconstructed pertur- 
batively (see eq. (0) below). f(x,l ± ) denotes the unin- 
tegrated gluon density, and x = (Q 2 + M 2 Q )/s and l± 
are the longitudinal momentum fraction and the trans- 
verse momentum, respectively, which are carried by a 
gluon inside the nucleon ^ ||. We substitute eqs. (|^) 
and (|J) into eq. (^), and go over to the "6-space" con- 
jugate to the fc_L-space via the Fourier transformation 
(j)(a,k±) — Air J d 2 be~ l ^- L ' b (f)(a,b\j b denotes the trans- 
verse separation between Q and Q. Then, it is straight- 
forward to see InxM ( - ±1 ' ) = which reflects helicity con- 
servation in the high energy limit. ImAl' -* can be calcu- 
lated in parallel with previous works for the vector meson 
production p|, [s[ pj|. In the integrand, there appear the 
terms e d± b and e~ d± b . The leading ln(Q 2 /A| CD ) con- 
tribution comes from the region Aq CD <C l\ <C Q 2 g{= 
[Q 2 +M 2 Q }/4) of the ^-integral §, §, where l^ b < 1 is 
satisfied because \b\ ~ 1/mg. We retain only the leading 
nonzero term in the power series in lj_ • b, which corre- 
sponds to the "color-dipole picture" jjj, [t9). Then, using 

J Qlt! dllf(x,ll)/ll = xG(x,Q 2 eff ) with G(x,Q 2 
conventional gluon distribution, we get (b = |6|) ui 



the 



. M{0)= -V2^W\ jwmQCA Qs(Qe2ff) 



xxG(x,Q 2 ef{ )[-£ 
Jo "(1 



N c M VQ Q cos aw 
1 daQ. 







2 dlnx 



a) 



dbb'PiaSKiibQm), (7) 



where Q m = [a(l — a)Q 2 + uiq] 1 ^ 2 , K\ is a modified 
Bessel function, and we have included the real part of 
the amplitude as a perturbation |5[ ||, ||. As expected, 
the result (0) is proportional to ca, so that the t\q meson 
is generated by the axial- vector part of the weak current. 

For comparison, we also calculate the diffractive vector 
meson production via the weak neutral current, by the re- 
placement ^l Q x , -> in eq. (§). Substituting eq. (Q), 
we find that £ = vj, i.e., Ai^ and M.^ 1 ^ give the pro- 
duction of the longitudinally and transversely polarized 



vector mesons, respectively, and that all these amplitudes 
are proportional to the vector coupling cy . (When we re- 
place the factor gwev / (2 cos 9w) by e c — 2e/3 in these 
results, we obtain the amplitudes A4^ identical to those 
for the diffractive electroproduction of J/ip ji], [j| , up to 
the corrections due to TZ of eq. 

Combining eq. ([?]) with the Z boson propagator and 
the weak neutral current by a neutrino, the forward dif- 
ferential cross section for the r\Q production is given by 



d 3 a{vN -> v'N't)q) 



dsdQ 2 dt 



9w 



t=o 



a{&kYeim 2 n s cos e w 



Q 2 



(Q 2 + M 2 ) 2 1 - e 
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where E v is the neutrino beam energy in the lab sys- 
tern, and e = [4(1 - y) - Q 2 /E 2 }/[2{1 + (1 - y) 2 } + 
Q 2 /E 2 ] with y = s/{2Mj^E v ) is the polarization pa- 
rameter of the virtual Z boson [jj], In order 
to evaluate the corresponding elastic t]q production 
rate, we assume the i-dependence as d 3 a/dsdQ 2 dt = 
d 3 a/dsdQ 2 dt\ t =o exp(Br, Q t) with a constant diffractive 
slope, as in the case of the vector meson production. 
Integrating d 2 a/dsdQ 2 = (l/B TlQ )d 3 a/dsdQ 2 dt\t=o over 
Q 2 and s, we get the elastic production rate cr(vN — > 
v'N'rjo) as a function of E v . 

For numerical computation of a(yN — > v'N'rjn), we 
need explicit form of (f>(a, b) of eq. (m). From eqs. (0) and 
(||), we can use the corresponding nonperturbative part 
of the vector-meson WF which was constructed in the 
previous works based on, e.g., non-relativistic potential 
model for heavy quarkonia |4|, [l9| . In this work, we adopt 
the WF from the Cornell potential model with the cor- 
responding quark masses m c = 1.5 GeV, m& = 4.9 GeV 
pfj. A procedure to construct <p(a, fcjj with the light- 
cone variables from the usual Schrodinger WF <pNn(k) 
has been given in Refs. 0, [nj. Also, we use the empir- 
ical values for the masses M Vc = 2.98 GeV, Mn — 0.94 
GeV and M z = 91.2 GeV. For r) b , we use M, )b = 9.45 
GeV estimated in Ref. [pi[ . Because the slope B VQ in- 
troduced above is unknown, we assume that B VQ has the 
same value as that for the corresponding vector meson: 
B Vc = 4.5 GeV~ 2 from the experimental value for J/tp 
Q and B Vb = 3.9 GeV~ 2 following the value for T in 
Ref. @. For the gluon distribution G(x, Q 2 ^) of eq. (g), 
we employ GRV95 NLO parameterization p3[ |. 

We show the elastic r/ c production rate, a(vN — > 
u'N'rjc), by the solid curve in Fig. |^. The result mono- 
tonically increases as a function of the beam energy E v . 
Such behavior is similar with that observed in the 7r- 
production data by the neutral and charged currents 
||. For comparison, we show the elastic J/tp produc- 
tion rate a(vN — > i/N'Jhp) by the dashed curve (see 
the discussion below eq. (0)) p3|. The rate for t] c pro- 
duction is much larger than that for J/tp by a factor 
~ 20. This is mainly due to the relevant weak cou- 
plings, ca for rj c and cy for J/ip, as (c^/cy) 2 = 7. 
Another important effect comes from the different be- 
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FIG. 3: The elastic production rates as functions of E v . Solid, 
dashed and dotted curves are for r\ c , J/tp and 77b , respectively. 
Note that the dotted curve shows the rate multiplied by 10 2 . 

havior of the Z light-cone WFs (Q) between the axial- 
vector and vector channels: It is straigtforward to see 
that the rj c diffractive amplitude (^) as well as the cor- 
responding amplitudes for the longitudinally and trans- 
versely polarized J/ip decreases rather rapidly for high 
<2 2 region, Q 2 > Am 2 ., so that the elastic production 
rates are in principle dominated by the low Q 2 contri- 
bution of the differential cross sections (similar behavior 
is observed also for J/tp electroproductions j|, ||). F° r 
weak neutral current proceses, however, a typical factor 
Q 2 /{Q 2 + Mg) 2 appears in the diffrential cross sections 



as in eq. (||) . The low Q 2 contribution of the differential 
cross section for J/ip is suppressed by this factor, while 
that for r\ c avoids the suppression due to ~ 1/Q behavior 
of eq. (0). This particular behavior of the r\ c diffractive 
amplitude comes from the spinor matrix element of the 
Z light-cone WFs (|) as u x {k)^ lb v v {q - k) ~ 1/Q for 
Q -> 0. On the other hand, u x (k)^°h x ,(q - fc) ~ Q 
and u\(k)^ ±1 " > vxi (q — fc) ~ const, for the vector channels 
relevant for J/ip. 

In Fig. |[ we also show the rjb production rate a(vN — > 
v 1 N'r)b) by the dotted curve. Although the rate for rjb is 
generally much smaller than that for rj c , the former in- 
creases more rapidly than the latter for increasing E v . 
Therefore, the rjb production rate could become compa- 
rable with the rj c production for higher beam energy, sug- 
gesting a possibility to observe rjb through the diffractive 
productions by high-intensity neutrino beams available 
in ongoing or forthcoming neutrino facilities. 

In conclusion, we have computed the diffractive pro- 
duction cross sections of i] c j, mesons via the weak neutral 
current, using the new results of the light-cone WFs for Z 
and r\ c ,b- Our results demonstrate that neutrino- induced 
productions will open a new window to measure rj c jy. 
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